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Abstract 

The Parisi formula for the free energy in the spherical models with mixed even 
p-spin interactions was proven in Michel Talagrand [16] . In this paper we study the 
general mixed p-spin spherical models including p-spin interactions for odd p. We 
establish the Aizenman-Sims-Starr scheme and from this together with many well- 
known results and Dmitry Panchenko's recent proof on the Parisi ultrametricity 
conjecture [IT] , we prove the Parisi formula. 
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1 Introduction and main results. 

In the past decades, the Sherrington-Kirkpatrick model [13] of Ising spin glass has been 
intensively studied with the aim of understanding the strange magnetic behavior of certain 
alloys. Despite only being a mean-field model, the formula for its free energy has an 
intricate and beautiful variational formula, first discovered by G. Parisi |Hl [9] , based on a 
hierarchical replica symmetry breaking scheme. As it is hard to compute the free energy 
from the Parisi formula, the closely related spherical model was studied by A. Crisanti 
and H. J. Sommers [I]; this model is widely believed to retain the main features of the 
Ising case, and the analogous Parisi formula admits a more explicit representation for the 
free energy. 

Mathematically, the Parisi formula for the Ising and spherical spin glass models was 
proven rigorously and generalized to mixtures of p-spin interactions in the seminal works 
of M. Talagrand [T6l [T7] following the discovery of the replica symmetry breaking inter- 
polation scheme by F. Guerra [7]. For technical reasons, only the mixture of even p-spin 
interactions was considered. The main missing ingredient in proving the general case was 
the validity of the ultrametric property for the Gibbs measure, which was recently veri- 
fied by D. Panchenko [UJ. This together with a series of known results led to the Parisi 
formula for the Ising spin glass model with general mixed p-spin interactions 
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In view of the approach in [12] , the most crucial step is the derivation of the inequality 
that the lower limit of the free energy is bounded from below by the Parisi formula. It 
is well-known, in the case of mixed even p-spin models, that the limiting free energy can 
be represented by a new variational formula through the Aizenman-Sims-Starr (A.S.S.) 
scheme pQ. An important implication of this scheme is that it provides a natural lower 
bound for the lower limit of the free energy for general mixed p-spin models. Since under 
a small perturbation of the Hamiltonian the asymptotic Gibbs measures in the A.S.S. 
scheme satisfy the Ghirlanda-Guerra (G.G.) identities, the main result in [TT] implies 
that the support of these measures are asymptotically ultrametric. This property then 
allows us to characterize the asymptotic Gibbs measure via the Poisson-Dirichlet cascades 
and consequently recover the Parisi formula in the general mixed p-spin case. 

In the present paper, we will establish the Aizenman-Sims-Starr (A.S.S.) scheme and 
prove the Parisi formula in the spherical model with general mixed p-spin interactions 
taking the approach of [12]. However, the task turned out to be much more difficult, 
mostly, due to the fact that uniform measure on the sphere is not a product measure. 
In particular, as we shall see in the proof of Theorem below, the establishment of the 
A.S.S. scheme in this case becomes much more intricate since in general its construction 
requires a certain decoupling of the configuration space. Also, the approximation of the 
Parisi formula using the ultrametricity of the Gibbs measure will be more involved. These 
obstacles are the main concerns of the paper and will be studied in the next two sections. 

Let us now state the mixed p-spin spherical model. For given N > 1, we denote by 
Sn = {cr G M. N : ||er|| = the configuration space and by A at its normalized surface 

measure. Consider the pure p-spin Hamiltonians Hn, p ((t) for p > 1 indexed by cr e Sn, 

HnA* 7 ) = N ( P -i)/2 9h,...,i v ^h---^i v , (1-1) 

l<ii,...,i p <jV 

where the random variables gi lt ...,i p for all p > 1 and all {i\,...,i p ) are i.i.d. standard 
Gaussian. The Hamiltonian of the mixed p-spin spherical model is given by 

H n {*) = Y,P p Hn, p {<t) (1-2) 

p>i 

with coefficients ((3 P ) that satisfy (3 P > and decrease fast enough, J2 P >i ^ v Pp < 00 • This 
technical assumption ensures that on the series (11.2j) converges a.s. and the covariance 
of the Gaussian process is given by a function of the normalized scalar product, called the 
overlap, i?i j2 = N~ l J2i<N \ <j 1 °f s Pi n configurations cr 1 and cr 2 , 

^H n {ct 1 )H n {cj 2 ) = N£(R h2 ), 

where E is the expectation with respect to all Gaussian randomness and £(x) = ^2 p>1 (3 p x p . 

Now we formulate the Parisi formula. Given k > 1, let m = (m^)o<£<fc and q = 
(<2^)o<£<fc+i be two sequences satisfying 

= mo < mi < • • • < irik-i — m k = 1; 
= go < qi < ■ ■ ■ < Qk < Qk+i = 1. 
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One may think of the triplet (k, m, q) as a distribution function x on [0, 1] with x(q) = rri£ 
if 1i < 1 < Qe+i for < £ < k and x(l) = 1. We define for 1 < £ < fc, 

e<p<k 

For a given parameter 6 > di, we define D(_ = b — for 1 < £ < k and D^ + i = b. Let 
#(x) = — £(x). The formulation of the Parisi formula involves a functional V of the 

distribution function x, 

P(x)= inf 1(6-1 -logfe + ^+ E -^S^- E m ^+i)-^*)) 

which is the limit of the Parisi functional described below. For every M > 1, let z p — 
(z P} i, . . . , z p< m) for < p < k be centered Gaussian random vectors with Ez Pti z Pt i> = 
5i t i'(£,'(q p+ i) — (,'(q p )) for any i ^ %' and let (z p )o< P <fc be independent of each other. 
Starting from 

X^ +1 = log / exp e ■ (z H h z k ) d\ M (e), 

J Sm 



we define decreasingly for < p < k 

K = — log E p exp mpXp^ 



rri p 



where E p denotes the expectation in the random vectors (zi)i> p . When m p = 0, this means 

3 Parisi functional 

V M {x)=X^- \ E m p (9(q p+1 )-e(q p )). 



Xf = EpX^. Set the Parisi functional 
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i< P <fc 



A series of nontrivial applications of elementary large deviation principle implies V(x) = 
liniM-s.oo Vm(x) (see Proposition 3.1 [IS]). Let us denote the partition function of the 
model by 

Zn = / exp H N (a)d\ N (cr). 

J Sat 



Theorem 1 (The Parisi formula). We have 

lim ^-E\ogZ N = inf V(x), (1.3) 
where the infimum is taken over all possible choices of x. 

The quantity in the limit on the left-hand side of (II. 3p is usually called the free energy of 
this model and the formula on the right-hand side of (II. 3p is the famous Parisi formula in 
the spherical form, which is also known to have the Crisanti-Sommers representation [I] 
(see Theorem 4.1 in [IS])- Of course, one may also consider the model in the present of 
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external field and prove the equation (II. 3p . but for the clarity of our discussion, we will 
only consider the Hamiltonian of the form (jl.2p . 

Acknowledgement. The author would like to thank Dmitry Panchenko for motivating 
this work and his advisor, Michael Cranston, for many helpful discussions. Special thanks 
are due to Alexander Vandenberg-Rodes for several helpful suggestions regarding the 
presentation of the paper. 



2 The A.S.S. scheme for spherical models. 

The main result of this section is the construction of the A.S.S. scheme for spherical 
models. For M, N > 1, let us denote by (z(cr) = (zi(er), . . . , zm(&)) ■ & G Sn) an d 
(y(er) : cr G Sn) two independent families of centered Gaussian random variables with 
covariances 

Ez i {(T 1 )z i/ {(T 2 )=5 i>i ^{R 1>2 ) (2.1) 
Eyicr'M* 2 ) = 6(R 1>2 ) (2.2) 

for every cr 1 , cr 2 G Sn- These random variables are independent of all the other random- 
ness. Let Gm,n be the Gibbs measure corresponding to the Hamiltonian 

H m ,n{°) = X) (M + N)(p-W X 9ii,..,i P °h ■ ■ ■ Vi? (2.3) 

p>l V ' l<h,...,i p <N 

with reference measure A at on Sn and let (-) MN be the Gibbs average of Gm,n- The 
following is our main result. 

Theorem 2 (The A.S.S. scheme). For any 5 > and M > 1, the lower limit of 
X 1 7. 1( >g Z\ is bounded from below by 

-J-liminf ( Elog / (expe ■ z(a)) M N d\ M (e) - Elog (exp v / My(cr)\ ) 

M N-+oo \ J Sm ' \ / M,Nj , 2 A , 

where vm is the probability measure of the M -dimensional Gaussian density and 

A s = {e G R M : M < \\e\\ 2 < M(l + 5)}. 

Note that a standard large deviation result yields liniM^oo M -1 log Um(As) = for every 
5 > 0. Thus (12.41) gives a lower bound for the lower limit of the free energy that will 
play an essential role in proving the Parisi formula in Section |3j Suppose that //jy is 
the uniform probability measure on = {-1, +1}^. If we replace the measure space 
(Sn, Xn) by A*jv) an d ignore the error term —8£'(1) + M _1 log u M (As) in (12 ,4p . then 
Theorem [2] gives the A.S.S. scheme for the Ising spin glass model with general mixed 
p-spin interactions (see e.g. Section 15.8 [IS]). 
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We will take the spirit in [TJ to establish the A.S.S. scheme. However, since uniform 
measure on the spherical configuration space in the spherical models is not a product 
measure, we have to seek a proper way to approximate this measure by another one that 
possesses a certain product structure. To deal with this obstacle, we will use the idea 
of Poincare's lemma (see page 77 [H]) that the projection on any fixed M coordinates 
of the uniform probability measure Xm+n on the sphere Sm+n converge weakly to M- 
dimensional standard Gaussian measure as N tends to infinity. The argument of our 
proof proceeds as follows. Recall the definition of H m +n from (11. 2p . For p G Sm+n, let 
us write p = (<x, e) for the first iV coordinates cr = (ax, ... , a at) and the last M cavity 
coordinates e = (ex, . . . , £m) and write 



H m +n(p) = H m ,n((t) + ^ EiZi{a) + 7(p) 



i<M 



where H MN is given in (12.31) . the term SiZ^cr) consists of all terms in H m +n with only 
one factor £j from e present, and the last term is the sum of terms with at least two factors 
in e. From now on, without causing any ambiguity, for any given vectors u), uj 1 G R K for 
some K > 1, we define R(ui, uj') = K^ 1 X^x^ 6 ^- One can check easily that 



(2.5) 



EZ^Z^ct 2 ) = £' (j/^Ri* 1 , ^)) , (2.6) 

M 2 

E 7 ( P ) 2 < — _ — ( 2 - 7 ) 

Let us notice an elementary fact that for a real- valued sequence (cjv), if lim inf at^oo c^/N = 
rGl, then M -1 liminfjv-s.oo(cAf+7v — cn) < t for every M > 1. This implies 

lim inf — — > — lim inf E log M+N (2.8) 

N^oo N ~ M N^oo Z N K ' 

and thus to obtain Theorem [21 it suffices to prove that (12. 4p is the lower bound of the 
right-hand side of f 12 . 8 j) . Let us write Elog Zm+n/Zn as 



E log j (exp e- z{cr)) MN dv M (e) - E log y 



f SN expH N d\ N (cr) 



+ Elog 



ISn 

7j M+N 



expH M ,Nd\ N (cr) 



(2.9) 



J A JM,Ndl>M{£) 



where 



Jm,n(s) = / exp(H M , N ((r) + e ■ z(cr))dA 7V (cr), e G W Ai 



In what follows, we prove several lemmas that will be used to bound these three terms in 
(12. 9p as iV tends to infinity. 
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Lemma 1. For z E M M , we have 



/ expe ■ zdu M (s) > u M (A s ) / exp e ■ zd\ M {e). ( 2 -10) 

J A s JSm 

Proof. Using the rotational invariance of i>m, there exists a probability measure 7m on R + 
such that v M is the image of 7m x Am under the map (s, e) i-> se. For 1 < s < (1 + S) 1 ^ 2 , 
using symmetry, 

J exp se ■ zd\ M (e) = J cosh(s£ • z)d\ M (e) 

> J cosh(e ■ z)d\M{s) = J exp e ■ zd\M{s). 

Therefore, (12.101) holds by 



r Hi+<5) 1/2 r 

/ exp e ■ zdv M [e) = / / exp se ■ zd\ M (£)d^ M {.s) 

Ja s Ji Js M 

Hl+<5) 1/2 r 

> / expe • zdX M {£)d^ M {s) 

Ji Js M 

= 7m({1 < s < (1 + 5) 1/2 }) f expe- zdX M {e). 

J Sm 



From this lemma and Fubini's theorem, the first term in ( 12. 9j) is bounded from below by 



Elog / (expe ■ z(cr)) MN dX M (e) + \ogu M (A s ) 



(2.11) 



Using the standard Gaussian interpolation technique, Lemma |2] below takes care of the 
second term in ( 12. 9p . 

Lemma 2. We have 

L expH N (a)d\ N (cr 



lim 

N-*oc 



Elog 



r a 1 \j\ 1 \ -Elog(exp V My(a 

J s expH MjN (cr)dX N ((T) \ 1 .\/..\ 



(2.12) 



Proof. Suppose (y'(cr) : a e SV) is a family of centered Gaussian random variables 
independent of all the other randomness such that in distribution 



H N {a) = H M}N (<t) + VMy'(cT). 
It is easy to check that Ey' (a x )y' (a 2 ) = J2 P >i ^D p N p R(cr l , cr 2 ) p , where 



D 



1 



• " M \Np^ (M + N)p~ 1 
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satisfies N P D P —tp — 1. Recall y(<r) from (12. 2p . Let us consider the interpolating free 
energy 



tp{t) = E log (^exp \fM (yiy'{cr) + y/1 - ty{a) 



Then |^(1) — ^(0)| gives the quantity in the limit on the left-hand side of (I2.12p . Thus 
it suffices to prove that ip'(t) — > uniformly in t G (0, 1) as N — > oo, which can be easily 
verified by using Gaussian integration by parts, as N — > oo, 

sup \<4>'(t)\ < M sup lEy'icr^y'i* 2 ) - Ey^M* 2 )] 



M sup 



p>i 



-»■ 0. 



Let us continue to find a lower bound for the lower limit of the third term in (12. 9p . It relies 
on an elementary lemma that allows us to decouple the surface measure (Sm+n, ^m+n) 
as the product measure (Sn x A$,\ n x F m ) asymptotically For K > 1, we denote by 
Sk the unit sphere in WL K and by \Sk\ the area of S^. Set 



M 



A m,n = J] [~\/M + N + l-j, y/M + N+l-j 



J M,N 



M 

j=1 \ D M+N+l-j 



I I 
\°M+N-j\ 



i\y/M + N+l-j' 



Define 



Fm,n(£) 



F M (e) 



M , 

b M> N n ( 1 - 

1 v M/2 

) exp 



M + N+l-j 

11-112 



M-j-JV-j-2 
2 



£ £ A 



M,Ni 



e e 



Define ai = 1 and for 2 < £ < M + 1, 



at(£) 



IT 1 



1-e 2 



3=1 



M + N -j 



£ £ A 



M,N- 



Lemma 3. Suppose that f is a nonnegative function defined on Sm+n- Then we have 



f(p)dX 



M+N' 



M+N 



p) = F M)N {e)de / / (cra M+l , e&x, . . . , e M a M ) d\ N (er) 



(2.13) 
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Proof. This is an elementary calculus formula that follows by an induction argument 
starting from 



K-3 



f(p)d\ K (p) = -j^L ! * (l- £ l\ 2 de I f (W^4,£ 1 r/A/,-_i.(<r) 



'Sk \S 1 k \VKJ-Sk\ K) ) Sk _ x \ \ K-l 

for every K > 2. 

□ 

As an immediate consequence of Lemma [3], the cavity coordinates converge weakly to 
the M-dimensional standard Gaussian random vector as N tends to infinity. We write 
ElogZ M+N / J As J MjN dv M (e) as 

Zm+n I a s Fm.n J M,Nds 

Elog M+N + JA, ■ _ (2 U) 

-T M,N>J M,NU>£ J As ^M'JM,nO'S 

Then the first term of (I2.14p can be controlled by the following lemma. 
Lemma 4. We have 

liminfElog T ^M±E >-M££'(l). (2.15) 

iv-^-oo F M>N J M ,Nde 

Proof. Suppose that (H' M N (cr) : cr G Sjv) is an independent copy of (Hm,n(c) '■ & G £jv) 
and is independent of all the other randomness. We consider the following interpolating 
Hamiltonian 

tp{t) =Elog / F M , N (e)de [ exp^^p) + # 2>t (p) + H 3>t {p))d\ N {a) , 
Ja s Js n 

where for p = (cr, e) with cr G Sn and e e Aj, 



Hi, t {p) = VtH M)N (o-a M +i) + Vl - tH' MN (c 

#2,t(p) = Z~L £i (yi a i Z i( a ' a M+l) + Vl - 



cr 



From (I2.13p . obviously <p(l) < K\ogZ M+N when iV is large enough. On the other hand, 
using Gaussian integration by parts and replicas p 1 = (cr 1 , e l ) and p 2 = (cr 2 , e 2 ), it follows 
that 



^ 2 

i<3 

where (-) t is the Gibbs average corresponding to random weights Fm,n Gxp(H^ t + H 2) t + 
i?3,t) with reference measure A at x de on SV x As and 

U x {p\p 2 ) = EH M , N (er l aM + i{e l ))H M , N (er 2 a M+ i(e 2 )) -EH'^^H'^a 2 ) 

U 2 (p\p 2 ) = J2 £l i £ i (^(e 1 )«i(e 2 )E^(o- 1 a M+1 (e 1 ))Z i (cr 2 a M+1 (e 2 )) - E^cr 1 )^ 2 )) 

i<Af 

U 3 (p l ,p 2 ) = E7(cr 1 a A f + i(e 1 ),£:iai(£ 1 ), . . .,4 fl M( £l )) 

• 7(cr 2 a A f+i(s 2 ),e 2 ai(e 2 ), . . . , e 2 v/ a M (e 2 ))- 
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From (123]) . (ESK, and (I277j) . it is easy to see that for j = 2, 3, 

lim sup \Ujip 1 , p 2 )\ -)> 

and from ( 12.5ft and using mean value theorem, 



lim sup [Z7 1 (/> 1 , p 2 ) | 

N -*°° P \p^S N xA 6 

< lim £'(L)(M + N) sup 



N 



a M+1 (E l )a M+1 (e 2 ) - 1 



sup 

£ 1 ,E 2 eA { 



M 



le 1 !!^ lie 2 !! 2 



M + iV 



Thus we conclude that lim sup sup 0<t<1 |<p'(t)| < M8£'(l) and so (12.15ft follows by 

Zm+n 



lim inf E log 



N 



> limin%(l) - (p(0)) 



> - lim sup sup \<p'(t)\ > -M8g(l). 

N^oo 0<t<l 



Finally, let us deal with the second term of f !2.14j) below. 
Lemma 5. We have 

J As F M N J M N de 
lim inf E log — ^ ■ — — > 0. 

AM-oo J F M J MyN d£ 



(2.16) 



Proof. Using the inequality log(l — x) > — 2x for < x < 1/2, when N is sufficiently 
large, we have for every e E As, 

Fm,n( £ ) 



log 



F M (e 
log 



M 1 "M + N-j-2 / 
&M,iv - - log - + ^ log ( 1 

7=1 V 



M + iV + 1 - j 



> log b M ,N ~ — log — + 



3||e| 



2 "°27r M + iV 

Note that limjv^oo bM,N = l/(27r) A/ / 2 as can be seen from (I2.13P with f — 1. By Jensen's 
inequality, we obtain 

f Ag F M J M>N exp (log F M>N / F M ) de 



lim inf E log 



> lim inf E 



f A F M J M>N log F M ^ N / F M de 



• v — * f A6 F M J MtN d£ 



3M 



> lim inf I log 6m a? log — , 



9 



and this gives (I2.16p . 

□ 

Proof of Theorem [H From (12T81) . fl2ll . and (I2.14p . using Lemmas 121 HI and El the lower 
limit of iV^ElogZ^v is bounded from below by 




- #(1) 



and then f l2~TTj) yields f l2~4j) . 

□ 

3 Proof of the Parisi formula for spherical models. 

This section is devoted to proving Theorem [TJ Our argument is started with Guerra's 
replica symmetry breaking interpolation scheme in the spherical models (see Section 3 
in [16]). In the case when p-spin interactions for odd p > 3 are not present in (11.21) . 
this scheme implies that the upper limit of the free energy is bounded from above by 
m{ x V(x). The fact that this bound still holds in the general mixed p-spin case relies on 
Talagrand's positivity principle on the overlap (see [15] and Section 12.3 in [18]) that is 
deduced from the validity of the G.G. identities [6] derived by a perturbation trick on 
the Hamiltonian. This part of the argument is well-known and has been used in several 
places. For example, the readers are referred to Theorem 14.4.4 [18] for a great detailed 
discussion. 

Our main concern in the proof will be concentrated on obtaining the lower bound, 
that is, the lower limit of the free energy is bounded from below by mi x V(x), that will 
be established in several steps and is based on the A.S.S. scheme. Again due to the 
non-product structure of the uniform measure on the sphere, additional technical issues 
occur. For instance, in the Ising spin glass model, the argument for the lower bound only 
uses the A.S.S. scheme with a single cavity coordinate, that is M = 1 (see Theorem 3.5 
[2], Lemma 11 [10], or Proposition 2 [12]), but this is not the case in our model. Indeed, 
recalling the formulation of the Parisi formula in the spherical model, one may figure out 
immediately that it is necessary to deal with the A.S.S. scheme (12 .4p with all M > 1 and 
compute their lower limit in M. Under this circumstance, we find the fact that the Parisi 
functionals (Vm)m>i have the same Lipschitz constant with respect to the Li-norm as 
stated in Lemma [7] below is of great use. 

The A.S.S. scheme for the perturbed Hamiltonian. For each u = (u p ) p >i with 
1 < Up < 2, let us consider a perturbation Hamiltonian 

p>i 

where H' Np (cr) are independent copies of the p-spin Hamiltonian in (11.11) . Since the 
perturbation term is of a small order, replacing H N with H N + H^ rt in (ll.3p obviously 
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does not affect the limit and in addition, the following A.S.S. scheme also holds. Consider 
another perturbation Hamiltonian 

ert 1 U f N \ {p - 1)/2 

H P m H n {ct) = (M + N y /S 2? \M + N ) H 'nJ^)- 

Let Gjj N be the Gibbs measure corresponding to the Hamiltonian H MN + H^ r N with 
reference measure and (•) w « be the average with respect to the product Gibbs measure 



M.N ■ 



Proposition 1 (The A.S.S. scheme). Recall A$ and vu from Theorem® For every 5 > 
and M > 1, we have 

lim inf — E log Z N 

> -J- lim inf sup (E log / (exps ■ z(cr)) MN dX M (e) - Elog (exp VMy(cr)) J 
m Af^oo u y y s \ / M,JV/ 



•5m 

1 

M 

(3.1) 



^'(1) + — logz/ M (A 5 ), 



where the supremum is taken over all u = (u p ) p >i with 1 < u p < 2. 

Starting with the Hamiltonian Hm+n + ^m+jv an d using the fact that H^+n * s °f a 
small order, one may easily find that the same arguments in the proof of Theorem [2] also 
work in Proposition (TJ so we will leave this part of the proof to the readers. In what 
follows, for each M > 1, we will construct a sequence (u m,n )n>i for u M,N = {u p u ' N ) p >\ 
with 1 < ul^ ,N < 2 along which the G.G. identities for the Gibbs measure Gj IN and the 
A.S.S. scheme for the Hamiltonian Hm,n + HjTn no ^ simultaneously. Let (a e )e>i be an 
i.i.d. sample from N and let R M ' N = (R^£i N )i,e'>i be the normalized Gram matrix, 
or matrix of overlaps, of this sample. Let C be the collection of all (p,n,f) with p > 1, 
n > 2, and / being a monomial of (Rf^ N )e,£'< n - We define 

<&M,N(p,n, f) 

for (p,n,f) G C. Then the G.G. identities can be stated as follows (see Theorem 12.3.1 

USD- 



Proposition 2. For M > 1, we have 

lim E u $ A f,Jv(p,n,/) = (3.2) 

/or a// (p,n,f) G C, where we view u = (u p ) p >i as a sequence of i.i.d. uniform random 
variables on [1,2] and E u zs the expectation with respect to the randomness of u. 
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Since C is a countable set, we can enumerate it as {(pj,rij, /j)}->i • Let us define 
$m,n(u) = ^j®m,n(Pj, Uj, fj), 

A m ,n( u ) = Jtf (^ E1 og J (expe ■ z((T)) MN d\ M (e) - Elog (exp VMy(cr)^ ^ . 

By Chebyshev's inequality, W u (*&m,n < ??) > 1 — ^u®m,n/v for every 77 > 0. If we pick 
77 = 2E u $M,Af, it implies that for any M and N we can always find u M,N = (Up I ' N ) p >\ such 
that $m,n(u m ' n ) < 2E u $ MjAr . Note that using $ MjN (n,p, f) < 2 and the G.G. identities 
( 13. 2 p implies limAr_*oo E u $m,jv = 0. Consequently, we have 

lim $m,n{u M ' N ) = (3.3) 
and from (13. ip . liminf7v->oo iV _1 Elog.ZAr is bounded from below by 

liminf A m ,n(u M ' N ) ~ + ^logu M (A s ). (3.4) 



In other words, for each M > 1, both the G.G. identities and the A.S.S. scheme hold 

r M,N 



along (u m ' n )n>i- Let us redefine the Hamiltonian H^ r N and the Gibbs measure 
by fixing parameters u = u M ' N . 

Asymptotic Gibbs' measures. Next, we will find a sequence of asymptotic Gibbs' 
measures that satisfy the G.G. identities and represent the lower limit in (13.41) . For each 
M > 1, let us pick a subsequence (N n ) along which the limits in (I3.3P and (13 ,4p are 
achieved and the distribution R M,N under KGjf^ converges in the sense of convergence 
of finite dimensional distributions to the distribution of some array R M . To lighten the 
notation, we assume without loss of generality that the sequence (N n ) coincides with the 
natural numbers. Under KG]^^ , the array R M,N is weakly exchangeable, that is, 



7r(f),7r(^) j ~ y%e 

for any permutation tc of finitely many indices. This property will still be preserved in the 
limit such that R M is a weakly exchangable symmetric non-negative definite array. Such 
an array is called a Gram-de Finetti array and is known to have the Dovbysh-Sudakov 
representation [5]. 

Proposition 3 (The Dovbysh-Sudakov representation). If{Ri^>)igi>\ is a Gram-de Finetti 
array such that Rg g = 1, then there exists a random measure G on the unit ball of a sep- 
arable Hilbert space such that 

(R e// ) ee!>1 =(T e -T e ' + 5 it e(l-\\T e \\ 2 )) 



Ll'>\ 



where (r £ ) is an i.i.d. sample from G. 
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Obviously the array R M satisfies Rf* e = 1. Using Proposition El let us denote by Gm the 
random measure generating the array R M and by (r ) an i.i.d. sample from Gm- We 
also denote by (-) M the average with respect to Gm- From (I3.3P and an approximation 
argument, the measure Gm satisfies the G.G. identities, 

1 n 

£=2 

for all n > 2, all bounded measurable function / of (i?^/)^/< n , and measurable function 
if). These identities will be used to identify the asymptotic Gibbs measures Gm in next 
subsection. 

Let us continue to check that the limit of Am,n(u m ' n ) in N can be expressed in 
terms of Gm, which relies on the following observation. Suppose for the moment that 
R = (Ri 7 £t)i,e r >i is an arbitrary Gram-de Finetti array such that Re/ = 1. Let C be its 
distribution and G be any random measure generating R as Proposition [3j Suppose that 
z(t) = (zi(r))i<M and y{r) are two Gaussian processes on the unit ball of some Hilbert 
space generating R with covariances 

E^t 1 )^ 2 ) = S^'ir 1 ■ r 2 ), Et/CtXt 2 ) = 6{r l - r 2 ), 

let r] be a M-dimensional standard Gaussian random vector, and let rj be a standard 
Gaussian random variable. Denote by (■) the Gibbs average corresponding to G and 
define 

* M (£) =^Elog / E„(exp e . (z(r) + 77(^(1) - C' (\\t\\ 2 )) 1/2 )) d\ M (e) 

J Sm 

- ^ElogE, (exp VM(y(r) +^'(1) - £'(||t|| 2 )) 1/2 )) - 

Lemma 6. For each M > 1, £ — > ^m{£) is a well- defined continuous function with 
respect to the weak convergence of the distribution C. 

This is a generalization of Lemma 3 in [12] and can be verified exactly in the same way. 
Let C M,N and C M be the distributions generating R M,N and R M , respectively. Then 
Lemma [6] leads to 

lim A m ,n(u M ' N ) = lim *m(£ M > N ) = *m(C m ). (3.6) 

Identitying the asymptotic Gibbs' measures using ultrametricity. Since the 
asymptotic Gibbs measure Gm satisfies the G.G. identities f)3.5p . the main result in [TT] 
implies that the support of Gm is ultrametric with probability one, that is, 

E(/(it!f 2 >mm«3,< 3 ))> M = l. 

This implies that for any q, the inequality q < r e ■ T e ' defines an equivalent relation and 
therefore, the array (I(q < Ree>))i,i'>i is nonnegative definite since it is block-diagonal 
with blocks consisting of all elements equal to one. Consider a sequence of functions 
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(F fe ) fe >i on [0, 1] by denning H k (q) = j/k if j/(k + l)<q<(j + l)/(k + 1) for < j < k 
and H k (l) = 1. Then we see that = {H k {Rf e )) t ^>i is symmetric, weakly exchangable, 
and satisfies the G.G. identities. In addition, it is non-negative definite since it can be 
written as a convex combination of nonnegative definite arrays, 

j= i \ ^ / 

Thus by the Dovbysh-Sudakov representation, Rff can be generated by a sample from a 
random measure G k M on the unit ball of some Hilbert space. Since sup 0<(?<1 \H k (q) — q\ < 
1/k, (R^) k >i converges weakly to R M . From Lemma El we then obtain 

lim 9 U (C?) = *m(£ M ), (3.7) 

where Cff is the distribution for R^ 1 . On the other hand, from the Baffioni-Rosati theorem 
[3], an ultrametric measure, such as G M , that satisfies the G.G. identities and under which 
the overlaps Rjf take finitely many values {p/k)o< p < k can be identified as the Poisson- 
Dirichlet cascade with parameters 

m M,k = {m M,k )o ^ k with m M,k = E <J ((jR M )i 2 < p/k) ^ M ) 
q M,k = {q M,k )o ^ k+1 with q M,k = p/{ J, + 1} _ 

For details of this result, one may consult the proof of Theorem 15.3.6 in [18]. Let x^ 1 
be the distribution function corresponding to this triplet (k, m M,k , q M ' fc ). Using the same 
argument as Theorem 14.2.1 in [18], we obtain 

^ M (Ct I )=V M (xt I ). (3.8) 

Lemma 7. For any two triplets (k, m, q) and (k', m', q') associated with the distribution 
functions x and x' respectively, we have 

\V M (x) - V M (x')\ <^ [ " *{q)\dq. (3.9) 

It is well-known that the same Lipschitz property as (13. 9p for the Parisi functional in 
the Ising spin glass model was established by F. Guerra [Tj. In the spherical models, 
despite of the non-product structure on the configuration space, one may find the same 
approach as Theorem 14.11.2 in [18] can also be applied to obtain (13. 9ft . Let J 7 be the 
collection of all distribution functions x on [0,1]. Define a metric d on T by letting 
d(x, x') = Jq 1 \x(q) — x'(q) \ dq for every x, x' G T . Using the Lipschitz inequality (13. 9ft . 
we can extend the functional Vm continuously to J~ ' . 

Proof of Theorem [TJ. If we denote by x M the distribution of Rf* 2 , it follows that by 
using (JX7D and (EBJ, 

^ M {C M )=V M {x M ). (3.10) 
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From Helly's selection theorem, (x m )m>i contains a weakly convergent subsequence with 
limit x°°. For clarity, we simply use (x M ) to denote this subsequence. Using ( 13.91) and 
f l3~T0]) . if we set x™ = H k (x°°), then 

limsup \^ M (C M ) - V M (xf)\ = limsup \V M (x M ) - V M (xf)\ 



2 

r(i) 



< ( limsup d(x M , x°°) + d(x°°,xf i 



2 

for every k > 1. To sum up, from (13. 4p . (13. 6p . and this inequality, we conclude that 

liminf ±-E\ogZ N > limsup (^ M (£ M ) + ±-\ogv M (A s )) - 

> limsup V M {x?) - limsup \* M (£ M ) - V M {x?)\ - 6?{1) 

M ->oo M -s>oo 

> inf V {x) - ^p-d(x°°, x%) - 8?(1) 

for every k > 1 and 5 > and this completes our proof by letting 5 — >• and k — > oo. 

□ 
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